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摘 要 :以 90° 大 弯 径 比 挠 性 弯 管 为 研究 对 象 ,基于 离散 相模 型 和 General 冲 蚀 模型 ,利用 Fluent 数 
值 仿真 软件 ,开展 不 同 工 况 下 气相 进口 速度 、 夹 带 固体 颗粒 质量 流量 、 弯 管 肋 条 及 弯 径 比 等 影响 因 
) 素 对 坷 弯 管 内 气 固 两 相 流 冲 蚀 特 性 的 数值 模拟 研究 。 研 究 结 果 表 明 :; 弯 管 纵 切 面 气 固 两 相 流 场 速 度 
分 布 存 在 分 离 现 象 ; 弯 管 的 最 大 冲 蚀 速率 随 气 相 进 口 速 度 的 增长 呈 指 数 关 系 , 随 颗粒 质量 流量 的 增 
长 呈 线 性 关系 , 随 肋 条 数目 的 增长 先 减 小 后 略微 增 大 最 后 趋 于 稳定 ;大 弯 径 比 弯 管 Sa 
强度 小 且 肋 条 的 存在 使 弯 管 的 耐 冲 蚀 性 能 得 到 一 定 的 提升 ; 弯 管 冲 蚀 最 严重 区 域 发 生 在 弯 管 沿 流 
人 左右 夹 角 的 区 域 且 呈 现 出 V 型 冲 蚀 形 貌 ,这 种 现象 是 由 于 弯 管 De 
迹 决定 的 ;在 弯 管 沿 流 动 方向 60° ~80° 夹 角 范 围 内 出 现 沟 楷 状 冲 蚀 形 貌 , 它 是 由 于 颗粒 在 弯 管 内 的 
二 次 及 多 次 碰撞 。 
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Study on erosion characteristics of gas-solid two-phase flow 
in a flexible elbow with a large bending diameter ratio 
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Abstract :This paper is intended for the study of a right-angle flexible elbow with a large bending diameter 
ratio. A numerical simulation on the erosion characteristics of gas-solid two-phase flow in the elbow is con- 
ducted employing the Fluent simulation software package ,in which such factors as gas inlet velocity ,mass 
flow rate of entrained solid particles ,elbow ribs and bending diameter ratio under different working cases 
are considered based on the discrete phase model and the Ceneral erosion model. The results show that the 


velocity distribution of the gas-solid two-phase flow field in the longitudinal section of the elbow is separa- 
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ted. The maximum erosion rate of elbow varies exponentially with the increase of the gas inlet velocity ,and 


linearly with the increase of the particle mass flow rate. As the number of ribs embedded in the inner wall 


of elbow increases ,the erosion ratio firstly decreases and then slightly increases and finally stabilizes. The 


strength of the secondary flow in the large bend is smaller than that in the elbow ,and the erosion resistance 


of the elbow improves to a certain extent due to the existence of ribs. The severely eroded area in the elbow 


occurs in the area where the elbow deflects along the flow direction at an angle of about 26° and presents 


V-shaped erosion morphology,which is caused by the special particle movement trajectory at the elbow. 


The groove-like erosion morphology appears within the angle range of 60° — 80° along the flow direction ， 


which is due to the secondary and multiple collisions of particles in the elbow. 


Key words :erosion wear; gas-solid two-phase flow;particle trajectory;maximum erosion rate 


目前 管道 运输 是 我 国 油气 运输 的 主要 方式 ,在 
天 然 气 开采 过 程 中 很 少 有 单 相 天 然 气 存在 的 情况 ， 
一 般 都 会 夹带 着 粉尘 砂粒 等 固体 颗粒 ,这 些 固体 颗 
粒 瑟 管道 中 输送 的 气体 构成 气 固 两 相 流动 体系 ,而 
固体 颗粒 等 杂质 对 管道 和 输 运 设备 的 安全 性 有 很 大 
的 谷 害 , 冲 蚀 磨损 是 造成 材料 破坏 和 结构 失效 的 重 
枚 原因 "0 。 大 弯 径 比 挠 性 弯 管 ( 弯 径 比 为 弯 管 曲率 
六 径 与 管道 直径 之 比 ,一 般 弯 径 比 >2.5 称 之 为 大 
索 黎 比 弯 管 ) 是 天 然 气 长 距离 输 运 中 不 可 或 缺 的 连 
接 伟 , 若 输 运气 体 夹带 颗粒 流 经 弯 管 等 局 部 位 置 时 ， 
固体 颗粒 会 对 管道 壁面 产生 冲 蚀 磨损 ,因此 对 管道 
营 管 局 部 位 置 的 冲 蚀 磨损 研究 具有 重要 意义 。 目 前 
计算 流体 动力 学 (CFD) 已 被 广泛 使 用 于 冲 蚀 磨 损 的 
数 御 模拟 计算 中 ,EDWARDS 等 中 利用 Ahlert 冲 蚀 
模型 对 索 径 比分 别 为 1. 5 和 5 的 90* 弯 管 进行 冲 蚀 
模拟 研 究 , 由 于 颗粒 与 管 壁 激烈 的 碰撞 导致 在 弯 管 
外 和 独处 产生 椭圆 形 的 冲 蚀 磨损 形 貌 。CHEN 等 ”以 
弯 软 比 为 1.5 的 标准 弯 头 为 研究 对 象 ,研究 其 冲 蚀 
形 貌 以 及 固体 颗粒 的 运动 轨迹 ,结果 表明 弯 管 处 冲 
仲 形 貌 为 V 型 。JASHANPREET 等 中 研究 了 浆液 输 
送 回路 弯 管 内 固体 颗粒 冲 名 ,数值 模拟 结果 表明 冲 
乌 率 随 流速 . 粒 径 和 管 径 的 增 大 而 增 大 。 较 低 的 冲 
蚀 磨 损 发 生 在 弯 管 入口 附近 并 在 弯曲 曲率 的 60° 附 
近 达 到 最 大 值 。 为 了 改善 气 固 两 相 流 管道 的 抗 冲 蚀 
性 能 ,目前 学 者 们 已 研究 出 多 种 方法 ,例如 SONG 
等 中 研究 发 现 ,在 直 管 段 内 壁 布置 一 定数 量 的 肋 条 
可 缓解 两 相 流动 中 管道 壁面 产生 的 冲 蚀 磨损 。 林 建 
忠 等 中 通过 实验 以 及 数值 模拟 研究 发 现在 管道 壁面 
上 开 纵向 较 浅 的 沟 槽 可 以 有 效 地 减轻 固体 颗粒 物 对 
其 壁面 产生 的 冲 蚀 磨损 。 季 楚 凌 等 "通过 建立 横 模 、 
纵 槽 和 凹 坑 等 仿生 表面 加 厚 弯 管 模型 ,经 数值 模拟 研 
究 发 现 采用 仿生 表面 技术 可 以 显著 改善 弯 管 的 抗 冲 


蚀 性 能 ,其 中 横 槽 的 抗 冲 蚀 性 能 最 好 。 文 献 [5-7 ] 
是 通过 改变 壁面 形状 来 影响 固体 颗粒 的 运动 轨迹 
从 而 提高 管道 的 抗 冲 蚀 性 能 。 一 些 研究 人 员 "* ”还 
发 现 设计 参数 也 是 影响 管道 冲 蚀 磨损 的 重要 因素 ， 
例如 弯 径 比 、 弯 管 曲率 和 弯 管 角度 。WANG 等 研 
究 了 不 同 弯 管 半径 下 的 弯 管 冲 蚀 磨损 ,结果 表明 增 
大 弯 管 曲率 是 减 小 由 冲 蚀 造成 弯 管 整体 磨损 的 有 效 
途径 。 

目前 很 多 研究 中 物理 模型 都 是 根据 无 颖 弯 头 规 
格 标准 建立 , 弯 径 比 大 约 为 1.5 的 管道 ,针对 大 弯 径 
比 弯 管 的 研究 较 少 。 本 研究 以 90" 大 弯 径 比 挠 性 弯 
管 为 研究 对 象 , 通 过 数值 模拟 的 方法 研究 挠 性 弯 管 
夹带 固体 颗粒 的 冲 蚀 磨损 性 能 ,分 析 影 响 弯 管 冲 乌 
性 能 的 相关 因素 ,并 与 小 弯 径 比 以 及 光滑 弯 管 进行 
对 比 , 探 讨 大 弯 径 比 以 及 挠 性 弯 管 的 特殊 结构 对 弯 
管 冲 蚀 性 能 的 影响 。 
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1 理论 计算 模型 


1.1 离散 相 控制 方程 


基于 拉 格 朗 日 坐标 系 对 连续 气相 中 的 离散 相 ， 
即 固体 颗粒 进行 运动 追踪 ,对 颗粒 运动 方程 积分 从 
而 得 到 颗粒 的 运动 轨迹 。Fluent 中 的 离散 相模 型 
(DPM) 忽略 固体 颗粒 之 间 的 碰撞 ,忽略 固体 颗粒 的 
破碎 , 故 只 适用 于 固 相 体积 分 数 小 于 10% 的 气 固 两 
相 流 计算 "。 颗 粒 运动 的 控制 方程 由 牛顿 第 二 运 
动 定律 得 出 , 即 
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其 中 : 为 电力 系数 ;v 为 连续 相 流 体 的 速度 ,m/s; 
uw 为 离散 相 颗 粒 的 速度 ,m/s;p, 为 颗粒 密度 ,kg/m ; 
p; 为 流体 密度 ,kg/m ; VPs 为 固 相 压力 ; /为 动力 符 
度 ,Pa. s; d, 为 颗粒 直径 ,m; C， 为 单 颗粒 表 观 中 力 
系数 ;Re 为 相对 雷诺 数 。 

在 气 固 两 相 流 场 中 ,颗粒 在 运动 过 程 中 不 免 要 
与 固体 壁面 发 生 接触 行为 ,这 必然 会 影响 流 场 中 颗 
粒 的 运动 轨迹 。 若 将 颗粒 认为 是 刚性 球体 ,与 壁面 
发 生 接 触 行为 时 要 在 接触 点 发 生 回 弹 ,其 回 弹 模型 
示意 图 如 图 1 所 示 。 颗 粒 以 人 射 速度 ww 向 壁面 运 
动 ,壁面 与 运动 方向 的 夹 角 9 称 为 冲击 角 , 颗 粒 在 接 
触 点 与 壁面 发 生 碰撞 出 现 能 量 损失 , 回 弹 速度 为 
v ,运动 过 程 中 能 量 的 损失 与 转移 使 得 颗粒 的 回 弹 
速度 低 于 入射 速度 。 恢 复 系数 。 用 来 表征 这 种 碰撞 
特性 , 即 碰撞 前 后 颗粒 速度 分 量 的 比值 ”。 本 研究 
选择 FORDER'" 的 法 向 恢复 系数 e, ,和 切 向 恢复 系 
来 描述 颗粒 与 壁面 的 碰撞 过 程 ,其 方程 如 下 


0.988 — 0.780 +0.190 -0.0240 +0.0270 
(3) 


Vu TV 


图 1 颗粒 与 壁面 碰撞 回 弹 模型 示意 图 


Fig.1 Schematic diagram of particle and wall 


collision rebound model 


1.2 冲 蚀 模型 


HUSER 等 研究 了 管道 中 的 砂粒 冲 蚀 磨损 ， 
并 提出 壁面 冲击 角 函 数 进行 冲 蚀 计算。 颗粒 冲 蚀 速 
率 可 以 在 壁面 边界 上 进行 监测 ,侵蚀 速率 定义 为 
VY mCCd Vf” 


- py 4 


p=1 


(5) 
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其 中 :7 为 壁面 冲 蚀 速率 ,kg .ms ;NN 为 碰撞 
颗粒 数目 ;m, 为 颗粒 质量 ;C(d,) 为 颗粒 直径 的 函 
数 ;7 为 颗粒 路 径 与 壁面 的 冲击 角 ;Ay ) 为 冲击 角 的 
函数 ;*" 为 颗粒 相对 速度 ;b(v) 为 颗粒 相对 速度 的 也 
数 ;4;.. 为 壁面 网 格 面 积 。Fluent 中 fy) 和 冲击 角度 
y 之 间 的 数值 关系 利用 线性 分 段 函数 描述 ,可 通过 
激 波 脉 冲 试验 获得 , 当 冲 击 角 yy 为 0°、20°、30°、 
45° .90° 时 ,对 应 的 冲击 角 函 数 A(y) 分 别 为 0、0. 8、 
1.0.0.5.0.4。 


2 数值 模拟 


2.1 几何 模型 建立 


含 砂 管 道 冲 蚀 实验 系统 由 离心 风机 、 砂 饶 .流量 
计 、 管 道 、 砂 箱 等 组 成 ,如 图 2 所 示 。 管 道 冲 蚀 磨损 
主要 发 生 在 挠 性 弯 管 部 分 , 故 本 研究 选取 90" 大 弯 径 
比 弯 管 作为 研究 对 象 进行 建 模 。 搁 性 弯 管 结构 如 图 
3 所 示 ,其 结构 由 互 锁 的 “S "型钢 带 组 成 , 挠 性 管 结 
构 具 有 良好 的 抗 坦 塌 性 与 耐 冲 蚀 性 。 


全) 和 砂 
亿 挠 性 过 管 
@ 砂 箱 


Fig.2 Schematic diagram of experimental system 
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到 3” 挠 性 弯 管 结构 示意 图 

Fig.3 Schematic diagram of flexible elbow structure 

根据 搁 性 弯 管 结构 建立 弯 管 内 流体 域 模型 进行 
数值 模拟 ,采用 肋 条 结构 来 近似 表征 流体 域 结构 , 几 
何 模型 由 人 口 直 管 段 六 、 弯 管 部 分 和 出 口 直 管 段 L， 
组 成 ,其 流体 域 模型 几何 尺寸 如 图 4 所 示 , 弯 管 进出 
口 管 径 D =25 mm, 弯 管 曲率 半径 尺 =200 mm, 弯 径 
比 RAD =8, 为 了 保证 流 场 充分 发 展 以 及 计算 结果 的 
准确 性 , 取 弯 管 进 出 口 直 管 段 长 度 为 6 倍 管 径 即 
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Li = 慷 =6D =150 mm。 肋 条 结构 几何 尺寸 如 图 5 所 
示 , 肋 条 截面 直径 D, =28 mm, 肋 条 宽度 h=2 mm。 
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弯 管 内 流体 域 几何 模型 图 


L g.4 Geometry model diagram of fluid domain in elbow 
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图 5 肋 条 结构 几何 图 


Fig. 3 Geometry of rib structure 
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边界 条 件 


”本 模拟 采用 Coupled 算法 ,流体 相 清流 模型 选 
用 Realizable 5-s 模型 ,利用 DPM 模型 模拟 颗粒 的 运 
动 轨迹 ,颗粒 相 湛 流 采用 随机 游 走 模 型 描述 ,采用 
Fluent 中 的 Ceneral 冲 蚀 模型 进行 冲 刨 计算。 
General 冲 刨 模型 的 参数 设置 如 下 :冲击 角 孔 数 
采用 分 段 函数 定义 , 当 冲 击 角 为 0* 、20° 、30°、45°、 
90° 时 函数 值 分 别 为 0、0.8、1.0.、0.5.0.4, 粒 子 直径 
函数 设置 为 常数 1. 168 4 x 10”, 速 度 指数 取 2.4。 
弯 管 入 口 边界 条 件 采用 速度 入 口 ,空气 作为 连续 相 
介质 , 进口 速度 大 小 分 别 取 为 10、15、20、25、30、 
35 m/s, 出 口 边界 采用 自由 出 流 , 壁 面 无 滑 移 。 离 散 
相 颗 粒 由 进口 面 释放 ,是 颗粒 的 初始 速度 与 流体 入 口 
速度 一 致 , 颗粒 相 在 壁面 处 的 边界 条 件 类 型 设 为 反 
弹 , 在 出 口 处 设置 为 逃逸 。 颗 粒 的 质量 流量 分 别 取 为 
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0.2.0.4.0.6.0.8.1.0 kg/s, 颗 粒 密度 为 2 650 kg/m ， 
颗粒 直径 为 0. 13 mm。 利 用 离散 相模 型 时 要 控制 颗 
粒 体积 分 数 小 于 10% , 当 颗 粒 质量 流量 为 1 kg/s， 
入 口 速度 为 10 m/s 时 ,颗粒 在 弯 管 中 的 体积 分 数 为 
7.69% , 故 所 有 工 况 满足 Fluent 离散 相模 型 计算 
条 件 。 


2.3 网 格 无 关 性 验证 


本 模拟 采用 非 结 构 化 网 格 , 为 了 准确 获得 弯 管 
壁面 及 肋 条 结构 处 的 冲 蚀 分 布 情况 ,对 肋 条 部 分 进 
行 局 部 加 密 。 为 了 得 到 合适 的 网 格 密度 保证 数值 模 
拟 的 精度 进行 网 格 无 关 性 验证 ,图 6 所 示 为 最 大 冲 
蚀 速 率 E,,, 与 网 格 数量 之 间 的 关系 曲线 图 。 由 图 
可 得 在 网 格 数量 较 小 时 , 弯 管 最 大 冲 蚀 速率 随 网 格 
数量 的 增 大 呈现 波浪 型 不 规则 变化 , 当 网 格 数量 达 
到 200 万 左右 后 , 弯 管 最 大 冲 蚀 速率 趋 于 稳定 。 
此 在 满足 计算 精度 的 条 件 下 ,尽量 减 小 模拟 计算 量 ， 
本 模拟 的 网 格 数量 为 226 万 ,网 格 质量 可 满足 计算 


4.0 
0 30 60 90 120 150 180 210 240 270 300 
n/10" 


6 最 大 冲 蚀 速率 随 网 格 数量 变化 曲线 


Fig.6 Curve of maximum erosion rate with the number of grids 


器 


2.4 模型 验证 


为 了 验证 模型 准确 性 ,与 VIEIRA 等 ”的 模拟 
结果 进行 了 对 比 ,如 图 7 所 示 。 以 文献 [17] 中 的 边 
界 条 件 为 依据 ,利用 Fluent 进行 气流 中 沙 粒 对 弯 管 
的 冲 蚀 模拟 ,气相 速度 取 为 23 m/s ,颗粒 直径 分 别 取 
为 150、300 hm, 膏 管材 料 为 316, 设置 密度 为 
7 980 kg/m 。 根 据 模拟 结果 弯 管 的 冲 蚀 形 貌 与 文献 
中 基本 一 致 ,最 大 冲 蚀 速率 与 文献 中 也 比较 接近 , 误 
差分 别 为 2.25% 和 6.5% ,在 误差 允许 范围 内 , 故 本 
研究 模型 可 行 。 


投稿 网 站 :http://cjam. xjtu. edu. en ” 微 信 公众 号 :应 用 力学 学 报 


二 应 用 力学 学 报 ChinaXiv 合 作 期 刊 第 40 卷 
Bh De 低 ,从 而 外 侧 流速 变 小 。 流 体 经 过 弯 管 内 测 时 由 于 
70le0 ws。 流体 的 惯性 边界 层 脱离 壁面 使 得 阻力 变 小 ,外 侧 流 
4.68e-05 2.43e-05 场 的 扰动 使 内 侧 压力 变 大 ,从 而 内 侧 流 速 变 大 。 截 
2342.05 人 取 弯 管 部 分 90* 截 面 上 的 速度 分 布 ,对 比 不 同人 口 速 
ooocrog 度 下 的 流 场 速度 分 布 ,如 图 9 所 示 , 随 着 流体 人 口 速 
[ees em) | 度 的 增 大 ,速度 分 离 区 域 有 变 大 的 趋势。 

300 hm 150 hm 


(a) 实际 模拟 结 


Erosion 
ey 


Erosion 
[kg .sim 


(b) 文献 模拟 结 
图 7 模拟 结果 对 比 图 


Fig.7 Comparison chart of simulation results 
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3 二 数值 模拟 结 果 与 分 析 


CD 

;内 连续 相 流 场 分 析 
ON 图 8 所 示 为 气相 进口 速度 在 10 ~35 m/s, 颗 粒 

直 移 为 0. 15 mm, 颗 粒 质量 流量 为 1 kg/s 的 条 件 下 

冤 生 纵 切 面 在 同一 标尺 下 的 速度 分 布 去 图 。 


inlet 
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弯 管 90。 截面 


图 8 不 同 入口 速度 下 弯 管 纵 切 面 速度 流 场 分 布 图 


Fig.8 Distribution of velocity flow field in longitudinal section 


of elbow at different inletvelocity 
连续 相 流 体 以 入 口 速度 稳定 地 流 过 弯 管 入口 直 
管 水 平 段 ,继续 流动 由 于 流向 发 生 改 变速 度 分 布 开 
台 出 现 明 显 的 分 离 现象 , 讨 管 外 侧 流 速 变 小 ,这 管 内 
侧 流速 变 大 。 流 体 流 经 弯 管 时 由 于 受到 离心 力作 用 
与 管 壁 发 生 碰撞 与 摩擦 ,造成 能 量 损 失 引起 速度 降 


(al) 10m 人 s (b)1SmMs (c)20m/s (d)25m/s (e)30m/s ({)35 m/s 
0.00 6.77 13.54 20.31 27.08 33.85 40.62 47.38 54.15 60.92 67.69 74.46 81.23 


Velocity m/s 
图 9 不 同 入 口 速度 下 弯 管 出 口 截 面 速度 云图 


Fig.9 Velocity cloud diagram of elbow exit section 


at different inlet speeds 
3.2 冲 蚀 影响 因素 


3.2.1 气相 进口 速度 

连续 气相 进口 速度 对 气 固 两 相 流 冲 蚀 速率 的 影 
响 如 图 10 所 示 , 对 应 工 况 为 和 人口 速度 从 10 m/s 逐 
渐 增 加 至 35 m/s ,颗粒 直径 为 0. 15 mm ,颗粒 质量 流 
量 为 1 kg/s。 从 图 中 可 看 出 弯 管 的 最 大 冲 蚀 速率 
E,, 随 着 进口 速度 "的 增 大 而 增 大 ,其 规律 近似 呈 指 
数 关系 ,说 明 冲 蚀 速率 主要 由 进口 速度 决定 。 当 气 
固 两 相 流 在 管道 中 流动 时 ,气相 流速 增加 ,由 于 颗粒 
具有 的 跟随 性 ,导致 回 相 颗粒 获得 的 动能 更 大 ,因此 
弯 管 的 最 大 冲 蚀 速率 变 大 。 
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图 10 气相 进口 速度 对 冲 刨 速 率 的 影响 

Fig. 10 The influence of inlet velocity on erosion rate 

图 11 所 示 为 弯 管 在 不 同 角 度 处 的 冲 蚀 速率 变 
化 图 , 沿 着 弯 管 流动 方向 分 别 取 0°、15°、30°、45°、 
60° .75° .90° 不 同 角 度 的 截面 进行 分 析 。 由 图 可 得 
弯 管 最 大 冲 蚀 速率 的 发 生 位 置 为 弯 管 沿 流动 方向 偏 
转 25° ~30° 区 域 范 围 内 , 冲 蚀 速率 出 现 第 一 个 峰值 ， 
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在 弯 管 沿 流动 方向 偏转 60° 左 右 达 到 二 次 峰值 ,此 时 
峰值 小 于 第 一 次 峰值 ,造成 这 种 现象 的 主要 原因 是 
由 于 颗粒 轨迹 的 改变 ,使 得 在 弯 管 的 不 同位 置 出 现 
不 同 程度 的 冲 蚀 磨损 。 对 于 同一 弯 管 角度 , 弯 管 的 
冲 蚀 速 率 随 着 和 人口 速 度 的 增 大 而 增 大 。 对 同一 人 口 
速度 , 弯 管 的 冲 蚀 速率 随 着 弯 管 角度 的 增加 有 一 定 
的 波动 ,出 现 峰值 ,但 总 体 上 呈现 出 减 小 的 趋势 ,这 
是 由 于 颗粒 与 讨 管 壁面 第 一 次 进行 碰撞 后 部 分 动量 
因此 损失 ,在 弯 管 特殊 结构 的 作用 下 颗粒 继续 与 壁 
面 发 生 二 次 及 多 次 碰撞 ,使 得 颗粒 对 壁面 的 冲击 变 
小 ,从 而 冲 蚀 速 率 变 小 。 


0/0) 
弯 管 不 同 角度 处 最 大 冲 蚀 速率 (不同 人口 速 度 ) 


.11 Maximum erosion rate at different angles of elbow 


( different inlet velocity) 

3 颗粒 质量 流量 

:= 颗粒 质量 流量 对 冲 蚀 速率 的 影响 如 图 12 所 示 ， 
里 ;质量 流量 从 0. 2 kg/s 逐渐 增加 至 1.0 kg/s, 气 
相 腓 口 速度 为 25 my/s ,颗粒 直径 为 0. 15 mm。 由 图 
可 得 弯 管 最 大 冲 人 名 速率 随 着 颗粒 质量 流量 m 的 
前 如 而 增 大 ,二 者 近似 为 线性 正 相关 。 颗 粒 质量 流 
量 的 增加 导致 相同 时 间 内 与 管道 发 生 碰撞 的 颗粒 数 
量 增多 ,因此 冲 蚀 速率 变 大 。 
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图 12 ”质量 流量 对 冲 蚀 速率 的 影响 


Fig.12 The influence of mass flow on erosion rate 


图 13 所 示 为 弯 管 不 同 角度 处 最 大 冲 蚀 速率 变 
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化 图 ,呈现 出 与 图 11 相似 的 规律 ,最 大 冲 蚀 速率 的 
发 生 位 置 为 弯 管 沿 流动 方向 偏转 25° ~30° 的 区 域 范 
围 内 ,图 中 出 现 的 两 个 峰值 对 应 弯 管 上 两 个 不 同 程 
度 的 冲 蚀 发 生 区 域 。 对 于 同一 弯 管 角度 ,在 弯 管 角 
度 45° 之 前 冲 蚀 速率 随 着 颗粒 质量 流量 的 增加 而 增 
大 ,但 在 弯 管 角度 45° 之 后 存在 冲 蚀 速率 随 着 颗粒 质 
量 流量 的 增加 而 减 小 的 情况 ,ZHENG 等 呈 研究 发 
现在 颗粒 质量 流量 增 大 至 一 定 程度 后 , 冲 蚀 速率 的 
增长 率 随 颗粒 质量 流量 的 增 大 而 减 小 ,因为 质量 流 
量 较 大 时 ,颗粒 间 碰 撞 的 概率 增 大 ,并 如 前 所 述 , 颗 
粒 速度 增 大 导致 碰撞 损失 的 动量 更 大 ,因此 与 弯 管 
壁面 碰撞 时 的 动量 减 小 从 而 使 得 冲 蚀 速率 有 所 降 
低 ,从 而 存在 同一 弯 管 角度 下 随 着 颗粒 质量 流量 增 
大 最 大 冲 蚀 速率 反而 下 降 的 情况 。 
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质量 流量 全 
一 是 一 2 TOON 
20| -04kgs 
二 ——0.6 kg/s | \ 
EE —s—0.8kgs | 人 
5 lO0kgs// 1 
‘mn 1 AN\ 1 
各 
Ed 
呈 
| 
由 5 


sl5 “0O 20 nn 7505 
0/(°) 


图 13 弯 管 不 同 角 度 处 最 大 冲 蚀 速率 (不 同 质量 流量 ) 


Fig.13 Maximum erosion rate at different angles of elbow 


(different mass flow) 
3.2.3 ” 肋 条 及 弯 径 比 
图 14 所 示 为 不 同 肋 条 数目 对 弯 管 最 大 冲 创 速率 
的 影响 , 弯 管 肋 条 数目 逐渐 从 26 增加 至 56 个 ,分 别 模 
拟 气相 进口 速度 为 20.25 m/s 时 的 弯 管 冲 蚀 情况 。 
4.5 
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图 14 弯 管 肋 条 数目 对 冲 包 速率 的 影响 
Fig.14 The influence of the number of 


elbow ribs on erosion rate 
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由 图 可 得 弯 管 最 大 冲 蚀 速 率 ,。 随 着 肋 条 数目 
的 增加 先 减 小 后 略微 增 大 最 后 趋 于 稳定 ,适当 增加 
肋 条 数目 可 以 有 效 减 小 弯 管 的 冲 蚀 速率 。 随 着 弯 管 
肋 条 数目 的 增加 ,管道 内 颗粒 与 壁面 碰撞 的 次 数 增 
加 ,颗粒 的 冲击 动能 被 进一步 减弱 ,从 而 使 得 冲 蚀 束 
率 下 降 。 随 着 肋 条 数目 的 继续 增 大 ,单位 时 间 内 颗 
粒 与 壁面 磁 撞 的 频率 增 大 , 冲 蚀 速率 略微 增 大 , 当 肋 
条 数目 增加 到 一 定 程度 后 发 现 再 增加 肋 条 对 弯 管 冲 
蚀 速率 的 影响 较 小 , 冲 蚀 速率 最 终 趋 于 稳定 。 

从 弯 管 冲 蚀 形 貌 可 以 得 出 冲 蚀 主要 发 生 在 弯 管 
外 侧 ,为 了 研究 其 具体 分 布 规律 取 弯 管 外 侧 壁 面 与 
中 心平 面 的 交 线 即 弯 管 外 侧 中 心 轮廓 线 上 的 冲 蚀 速 
率 ,如 图 15 所 示 得 出 弯 管 外 侧 轮廓 线 上 的 冲 蚀 速率 
沿 着 弯 管 流动 方向 偏转 角度 的 变化 情况 ,同时 对 比 
了 同 尺 寸 光 滑 弯 管 外 侧 中 心 轮廓 线 上 的 冲 蚀 速 率 分 


源 - 而 挠 性 雹 管 沿 着 外 侧 轮 廊 线 在 一 定 的 范围 内 上 
人 要 动 有 部 分 位 置 冲 他 速率 接近 于 0。 
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15 弯 管 外 侧 中 心 轮廓 线 上 的 冲 蚀 速 率 分 布 


Fig.15 Erosion rate distribution on the outer 


center contour of the elbow 

从 图 15 可 得 弯 管 最 大 冲 蚀 速率 发 生 在 偏转 角 
度 为 26 左右, 提取 弯 管 治 流动 方向 偏转 26" 附 近 一 
个 完整 肋 条 轮廓 线 上 的 冲 蚀 速率 分 布 如 图 16 所 示 。 
图 例 中 云图 为 弯 管 中 心平 面 上 的 冲 刨 速率 分 布 云 
图 ,从 图 中 可 看 出 肋 条 的 上 游 区 域 几乎 不 存在 冲 蚀 ， 
而 下 游 区 域 存在 明显 的 冲 蚀 区 域 , 这 与 弯 管 沿 肋 条 
轮廓 冲 蚀 速率 分 布 情 况 相 符合 , 冲 蚀 速率 在 肋 条 轮 
廊 线 上 先 增 大 后 减 小 ,最 大 冲 蚀 速率 发 生 在 肋 条 的 
下 游 区 域 , 即 肋 条 背面 与 弯 管 平 直 段 相 连 处 。 由 于 
肋 条 的 存在 挠 性 弯 管 壁面 凹凸 不 平 , 使 得 弯 管 上 存 
在 部 分 位 置 未 接触 到 流体 中 的 颗粒 从 而 使 得 冲 蚀 速 
率 为 0, 这 可 在 一 定 程度 上 组 减 弯 管 的 冲 蚀 。 
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图 16 弯 管 最 大 冲 蚀 速率 位 置 处 沿 肋 条 轮廓 冲 蚀 速率 分 布 


Fig. 16 Erosion rate distribution along rib contour at the 


position of maximum erosion rate of elbow 

取 进 出 口 管 径 D=25 mm ,曲率 半 径 R=37.5 mm， 
弯 径 比 RAD =1.5 的 小 弯 径 比 光滑 弯 管 在 同一 工 况 
下 进行 模拟 ,分 析 对 比 不 同 弯 径 比 弯 管 在 不 同 角 度 
处 弯 管 截面 流 场 速度 分 布 ,图 17 所 示 为 不 同 弯 径 比 
弯 管 沿 流动 方向 分 别 偏转 45°* 和 90° 时 弯 管 截面 速 
度 流 场 失 量 图 。 从 图 中 可 得 小 弯 径 比 光滑 弯 管 存在 
明显 的 二 次 流 现象 ,由 于 弯 管 中 心 流 体 速 度 较 大 从 
而 所 受 离心 力也 大 ,导致 二 次 流 以 较 大 的 速度 向 弯 
管 外 侧 流 动 ,管道 壁面 附近 的 流体 被 迫 向 内 侧 流 动 ， 
最 终 形 成 了 双 游 涡 二 次 流 分 布 。 相 比 于 小 弯 径 比 弯 
管 ,大 弯 径 比 弯 管内 部 流 场 二 次 流 现象 不 太 明 显 ,二 
次 流 强 度 较 小 。 
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图 17 不 同 弯 径 比 弯 管 在 不 同 角度 处 
弯 管 截面 速度 流 场 矢 量 氏 
Fig.17 Vector diagram of the velocity flow field of the elbow 


with different bending diameter ratios at different angles 
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从 图 17(b) 可 以 看 出 弯 管 壁面 附近 存在 低速 的 
小 洲 涡 ,这 是 由 于 肋 条 的 存在 使 得 弯 管 壁面 附近 的 
稳定 流 场 被 破坏 ,从 而 使 得 此 区 域 的 清流 强度 大 于 
光滑 壁面 ,改变 将 要 撞击 到 壁面 颗粒 的 运动 方向 使 
其 流入 中 心 流 场 区 域 ,同时 低速 的 洲 涡 流 场 会 对 颗 
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霍 志 远 ， 等 ; 大 弯 径 比 挠 性 守 管 内 气 固 两 相 流 冲 伺 特性 相册 inaXiv 合 作 期 刊 44 


由 图 可 得 不 同 速 度 下 的 弯 管 冲 蚀 形 貌 呈现 出 相 
同 的 规律 , 冲 蚀 磨损 主要 发 生 在 弯 管 外 侧 , 冲 刨 最 严 
重 位 置 在 弯 管 沿 流动 方向 偏转 26° 左 右 区 域 , 冲 蚀 磨 
损 呈 现 出 V 型 冲 蚀 形 貌 , 弯 管 处 颗粒 的 运动 轨迹 表 
现 为 :颗粒 与 弯 管 进行 碰撞 后 ,其 轨迹 分 为 两 条 ,从 


粒 的 冲击 起 到 一 定 的 缓冲 吸收 部 分 动能 ,减轻 了 颗 
粒 对 壁面 的 冲 蚀 ,从 而 流体 在 小 弯 径 比 光 滑 弯 管 中 
所 引起 的 流动 损失 更 严重 , 肋 条 的 存在 使 弯 管 的 耐 
冲 蚀 性 能 得 到 一 定 的 提升 。 


3.3” 冲 蚀 形 貌 及 颗粒 运动 轨迹 分 析 


图 18 所 示 为 气相 进口 速度 分 别 为 10 ~35 m/s， 
颗粒 直径 为 0.15 mm, 质量 流量 为 1 kg/s 时 的 弯 管 
冲 蚀 磨损 速率 分 布 云图 。 
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18 不 同 速度 下 弯 管 冲 蚀 速率 分 布 云图 


Fig.18 Cloud diagram of erosion rate distribution of elbows 


at different velocity 


而 出 现 V 型 种 蚀 形 得 ,这 与 PEREIRA SOLNORDAL 
等 "”” 发现 的 规律 一 致 。 在 V 型 冲 蚀 形 貌 之 后 出 
现 了 沟 槽 状 冲 蚀 形 貌 ,主要 发 生 在 弯 管 沿 流动 方向 
偏转 60° ~80° 夹 角 范 围 内 , 沟 槽 状 冲 蚀 形 貌 距 离 V 
型 冲 蚀 形 貌 较 远 , 且 冲 蚀 速 率 小 于 V 型 冲 刨 形 貌 。 

图 19 所 示 为 弯 管 中 颗粒 的 运动 轨迹 追踪 图 。 
由 于 弯 管 结构 的 特殊 性 ,在 运动 过 程 中 颗粒 会 与 搅 
性 弯 管 壁面 以 及 肋 条 发 生 多 次 碰撞 ,颗粒 的 速度 和 
运动 方向 都 会 发 生 改 变 , 故 颗 粒 在 弯 管 中 的 运动 规 
律 复杂 。 颗 娄 从 进口 进入 弯 管 后 ,如 图 19(b) 所 示 ， 
与 弯 管 壁面 发 生 直接 碰撞 ,产生 初次 碰撞 并 改变 颗 
粒 运 动 方向 ,从 图 19(a) 可 得 颗粒 与 壁面 发 生 接触 
碰撞 后 运动 方向 主要 治 着 两 个 方向 进行 ,这 也 与 VY 
型 冲 蚀 形 貌 相符 合 。 对 于 大 弯 径 比 弯 管 ,初次 碰撞 
是 造成 弯 管 冲 蚀 最 严重 位 置 的 主要 原因 , 故 V 型 冲 
蚀 形 貌 主要 是 由 于 颗粒 的 初次 碰撞 引起 。 经 过 初次 
碰撞 后 固体 颗粒 在 弯 管 内 反弹 形成 二 次 及 多 次 碰 
撞 , 这 也 是 沟 权 型 冲 蚀 形 貌 形成 的 原因 。 


(a) 三 维 轨迹 
匀 19 ”颗粒 运动 轨迹 追踪 图 
Fig. 19 Particle tracking diagram 


(b) 纵 面 轨迹 


4 结论 


本 研究 采用 离散 相模 型 和 Ceneral 冲 蚀 模型 ,对 
90° 大 弯 径 比 找 性 管内 气 固 两 相 流 的 冲 蚀 特性 开展 
了 数值 模拟 研究 ,主要 结论 如 下 。 

1) 大 弯 径 比 弯 管内 部 流 场 速度 分 布 存在 明显 的 
分 离 现象 , 弯 管 外 侧 流速 较 小 ,内侧 流 速 较 大 , 且 随 
着 人 口 速 度 的 增 大 速度 分 离 区 域 越 大 。 

2) 气 相 进口 速度 对 冲 蚀 速 率 的 影响 较 大 , 弯 管 
的 最 大 冲 蚀 速率 随 着 气相 进口 速度 的 增 大 而 增 大 且 
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二 者 近似 为 指数 关系 , 弯 管 的 最 大 冲 蚀 速率 随 着 颗 
六 质量 流量 的 增 大 而 增 大 且 二 者 近似 为 线性 关系 ， 
弯 管 最 大 冲 蚀 速率 的 发 生 位 置 大 约 为 沿 弯 管 流动 方 
向 偏转 26° 左 右 。 

3) 弯 管 的 最 大 冲 蚀 速率 随 着 肋 条 数目 的 增加 
先 减 小 后 略微 增 大 最 后 趋 于 稳定 ,适当 增加 肋 条 数 
目 可 有 效 减 小 弯 管 冲 包 速率 。 弯 管 外 侧 中心 轮 廊 线 
上 的 冲 蚀 速率 分 布 不 均匀 且 存 在 冲 蚀 接近 于 0 的 区 
域 ,这 是 由 于 肋 条 的 特殊 结构 导致 , 相 比 于 光滑 弯 
管 , 肋 条 可 在 一 定 程 度 上 组 减 弯 管 冲 蚀 。 

4) 相 比 于 小 弯 径 比 ,大 弯 径 比 弯 管 内 ,流动 的 二 
次 流 强 度 小 且 肋 条 的 存在 使 弯 管 的 耐 冲 蚀 性 能 得 到 
一 定 的 提升 。 


__5) 弯 管 的 冲 蚀 磨损 主要 发 生 在 弯 管 外 侧 , 冲 蚀 


= 
ey 


严重 区 域 发 生 在 索 管 26° 左 右 并 且 旦 现 V 型 促 蚀 
戏 。 在 弯 管 60° ~ 80* 范 围 内 出 现 沟 槽 状 冲 蚀 形 
二 时 冲 蚀 速率 小 于 V 型 冲 蚀 形 貌 。 
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